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We present a theory of the optical control of the spin of an electron in an InAs quantum dot. We 
show how two Raman-detuned laser pulses can be used to obtain arbitrary single-qubit rotations 
via the excitation of an intermediate trion state. Our theory takes into account a finite in-plane 
hole gr-factor and hole-mixing. We show that such rotations can be performed to high fidelities with 
pulses lasting a few tens of picoseconds. 

PACS numbers: 78.67.Hc, 03.67.Lx 



The spin of an electron in a quantum dot (QD) is cur- 
rently viewed as one of the leading contenders in the drive 
to develop a qubit from which a practical quantum com- 
puter can be constructed The optical manipulation 
of such spins is an area of tremendous contemporary in- 
terest and offers a route to the ultra-fast control of qubits 
required to implement quantum logic. 

In this paper we provide a theory of optically- 
controlled single-qubit rotations of a single electron spin 
in a self- assembled InAs QD. Such dots are at the fore- 
front of experimental effort in this field, and their promise 
for quantum computation purposes has been illustrated 
in a number of recent experiments. For example, the 
spin state of a QD electron has been prepared in a 
pure state to a very high-fidelity 0], and the exchange- 
interactions between optically-excited coupled dots have 
recently been mapped 

The theory we describe here shows that by using two 
pulsed, Raman-detuned lasers one can obtain arbitrary 
single-qubit rotations of the electron spin. The rotation 
proceeds through the virtual excitation of an intermedi- 
ate trion state. We show that, for realistic parameters, 
this can be accomplished with laser pulses lasting just a 
few tens of picoseconds. This is significantly faster than 
the decoherence time of the electron spin in such a dot, 
which has been measured at a few microseconds Q • Fur- 
thermore, since the trion state is only virtually excited, 
this technique avoids problems arising from spontaneous 
emission from this state. 

This work is an extension of that of Chen et al. [f| , who 
have presented a theory of single-qubit rotations for dots 
in which the heavy holes have a zero in-plane g-factor. 
Subsequent experiments Q have shown that, although 
this is the case in GaAs fluctuation dots Q , the ^-factor 
in self-assembled InAs dots is actually finite. This is im- 
portant because an in-plane magnetic field is vital to the 
rotation mechanism. Furthermore, in InAs dots we ex- 
pect significant hole-mixing, and this also needs to be 
taken into account. Finally, the analytic technique that 
we use here is somewhat different from that employed in 
Ref. [f| , and is capable of describing the qubit rotations 
to a much higher fidelity. 



FIG. 1: The four-level model of the electron-trion system 
in the Voigt basis consists of two Zeeman-split single-electron 
ground states |x±) with spins aligned in the x direction, and 
two trion levels |r±) with heavy-hole spins also in the x di- 
rection. The ground states are split by Zeeman energy E B 
and the trion states by E B . Arrows indicate allowed opti- 
cal transitions with H and V denoting two orthogonal linear 
polarizations. To obtain single-qubit rotations we pump tran- 
sitions Vi and Hi with a common Raman detuning of 6. 



I. SPIN-TRION QUBIT SYSTEM 

We consider a singly-charged self-assembled InAs QD 
with growth direction z. The spin of the electron trapped 
in the QD is our qubit degree of freedom and we will per- 
form rotations of this spin through the virtual excitation 
of an exciton within the dot. Figure [1] shows the four- 
level model that describes the pertinent features of this 
system. We apply a magnetic field in the x-direction. 
The Zeeman energy of a QD electron in this field is 
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where g% is the electronic 



g- factor, \xb is the Bohr magneton, B x is the magnitude 
of the field, and s% = ±1/2 corresponds to the electron 
spin. 

The heavy-hole component of the trion also splits un- 
der this field in InAs QDs [1, [1], and can be described 
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with a Zeeman Hamiltonian Tig 



where s x = ±1/2 are the eigenvalues of a pseudo-spin, 
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the components of which correspond to heavy-hole states 
aligned in the s-direction, and g x is the hole (/-factor. Re- 
cent measurements have given the magnitudes of these 
g-factors as \g%\ = 0.46 and \g x \ = 0.29 Q. For concrete- 
ness, we follow Ref. Q here and assume that both these 
g-factors are negative, but this is in no way essential. 
In our initial treatment we will neglect hole mixing and 
show later that it can be incorporated into the analysis 
through a set of only very minor modifications. 

The four levels of our model are then: the two elec- 
tron ground states with spins in the x-direction, \x±) = 
2~ 1 / 2 (I I) i | t)), where | J.) and | |) represent electron 
spins in the z direction; and the two trion levels, |t±) = 
2" 1 (I IT) - I Tl» (I 4) ± I ft)), where U> = | § , - § ) 
and | 1*1") = | | , | ) denote heavy-hole states also aligned 



in the z direction. 

Figure [T] shows the allowed optical transitions between 
these levels. Due to the splitting of the trion level, these 
transitions are linearly polarized. We have defined the 
polarization vectors in terms of cr± circular polarizations 
as V = 2- 1 /2( (T _ + <t+) and H = 2- 1 / 2 (< T _ - <t+) [|. 

To obtain our qubit rotations, we illuminate the system 
with two phase-locked laser pulses propagating in the z 
direction. We use one V-polarized pulse with frequency 
toy and time-dependent Rabi frequency f2y(t), and one 
H-polarized pulse with frequency ujh and Rabi frequency 
f2ff(£). With these lasers we pump the two transitions to 
the trion level |r+), labelled Vi and Hi in Fig. [IJ with a 
common Raman detuning of S. The model Hamiltonian 
in the basis \x— }, \r+), |r_)} is then 
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where E T is the trion energy, a is the relative phase of the A# = (g e x — g^)jiBB x = E C B + Eg. (2) 

two lasers and we have set h = 1. The time-dependence 
of fly and fin is understood. 

We set the frequencies of the two lasers as ujy = 
E T — ^£_b — S and u>u = E T + ^A^ — 5, where we have 
introduced 

Transforming to a rotating frame, we obtain the Hamil- 
= {g e x +g^)^ B B x = E%-E h B tonian 
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The terms oscillating with frequencies As and Tib de- 
scribe the off-resonant driving of transitions H2 and V2 
respectively, both of which involve the state |r_). In 
deriving our analytic solution for the behaviour of this 
system, it is these off-resonant transitions that we treat 
approximately and it is therefore the magnitude of the 
quantities and Tb that determine the accuracy of 
our description. 



II. ADIABATIC ELIMINATION OF TRION 
LEVELS 

We now seek an compact, approximate description of 
the evolution of the qubit sector and we shall proceed 
through the adiabatic elimination of the trion levels. This 
is appropriate for the situation we consider in which the 
detuning is sufficient that we only virtually occupy the 
trion level. The equations of motion for the wave function 
coefficients under the action of H. are 

ZCx-\- ' =: i^yG Cf-{- ^tj^G ^ Cf 
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The equations for c T ± can be rewritten as 
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We proceed by integrating by parts and assuming that 
the envelopes flv,H and thereby c x ± are slowly-varying 
functions of time. We thereby arrive at an approxi- 
mate expression for the trion amplitudes in terms of the 
ground-state coefficients: 
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The validity of these expressions are conditioned on the 
following adiabatic constraints 



\Qhc x -\ > 
\flnc x+ \ > 
\fl v c x 



(7) 



which mean that, as well as the laser amplitudes being 
slowly varying, the laser frequencies can not be too close 
to resonance with any the transitions. 

Assuming that these conditions are met, we can sub- 
stitute c T ± from Eq. ^ into the remaining equations of 
motion Q , and obtain a closed set of equations of motion 
for the qubit coefficient c x ±: 



dt 
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= h(t) 
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where h(t) is the effective Hamiltonian acting solely on 
the qubit space: 
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In writing this Hamiltonian, we have neglected further 
terms oscillating with frequency 2Ej 3 . These terms are 
expected to be negligible since, without them, the Rabi- 
frequency of the effective model is fly fljj/6, and thus 
provided that 
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these terms are rapidly oscillating in comparison and thus 
approximately self-average to zero. This effective two- 
level Hamiltonian h(t) forms the basis of our approach 
and provides a good account of the system as we will 
show. 



III. TIME-EVOLUTION OPERATOR 

Given that the qubit evolves under the action of effec- 
tive Hamiltonian h(t), we may once again make use of 
the adiabaticity of the system to obtain an approximate 
form for U(t), the time-evolution operator of the system. 
Firstly, let us specify that the envelopes fly (t) and fl B (t) 
have the same shape but different amplitudes and write 
flnit) = uCly(t), where v is constant in time. We can 
then write h{t) in the form 



h(t) = -X(t)n ■ a + pA. 



(11) 



and can drop the constant term fit. In general, both the 
forefactor A and the unit vector n depend on time, but 
since we are assuming the same shape for both pulses, 
the axis n remains fixed, and all the time dependence 
is contained in X(t). We can then approximate the time 
evolution operator as 



U{t) = exp{-i/2A(i) n ■ er} . 



(12) 



where A(t) = J_ dt'\(t'). The final output gate opera- 
tor is U = U(t — ► oo). 

The operator U represents an arbitrary single-qubit 
rotation. The rotation angle at time t is given by 
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where 



r 2 = S|(Ab-^) 2 + A|(<5 + S b )V 
+2{A B - S)(S + E B ) 

x [A B (S B + 2S)iy 2 - 26(6 + E B )] v 2 (14) 

Note that this expression for A(f) depends only on the 
"area-under-the-pulse-squared" , and not on the details of 
the shape of the pulse. This is typical of the adiabatic ap- 
proach we are pursuing here. For simplicity, we will con- 
sider a Gaussian pulse envelope, fly = ^4exp (— t 2 /2T 2 ), 
in which case the total angle is 
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The axis of rotation has components 
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n,3 = sin (3 



(15) 



(16) 



4 



0.9 



0.8 



0.7 



0.6 



0.5 




0.1 0.2 
<5(meV) 



I 

» 1 

I i 



0.3 



0.8 



. 6 



0.4 



0.2 




0.1 



0.2 



X 

\ \ 






* 


u 


1 / 


/I 

• \ 


1// 



0.3 



5(meV) 



FIG. 2: The fidelity JF, which reflects the degree to which 
our description approximates the exact evolution of the sys- 
tem, as a function of Raman detuning S. Here we consider 
a 7r-rotation about the 3 axis, corresponding to a change in 
relative phase between the two qubit levels, using pulses of 
durations T = 5 ps (green dash-dot), 15 ps (blue dash) and 
50 ps (black solid line). As is clear, for this type of rotation, 
excellent agreement (T close to unity) can be obtained with 
even very short pulses (T ~ 5 ps). The model parameters 
are g% = —0.46 and g% = —0.29 with a magnetic field of 8 T. 
Laser parameters were v — 0, with A chosen to give A = n. 



where 
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and a is the phase difference between the pulses. 

Let us now describe two specific rotations as illustra- 
tive examples. We consider first a rotation about the 3 
axis, which corresponds to a change in the relative phase 
of the qubit states. Only a single laser pulse is necessary 
and thus we set v = 0. The rotation is then precisely 
about the 3 axis, n = (0, 0, 1), by an angle of 
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A rotation with the axis in the 1-2 plane includes a 
transfer of population between the two qubit levels. This 
requires two lasers and, by setting the relative amplitude 
of the two pulses according to 



' Z B (A B -5) 
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we set n 3 = 0, and obtain a rotation axis exactly in the 
1-2 plane. The direction of this axis is then specified by 
the relative phase of the two lasers, a. In order for v to 
be real and finite we require that A B < 5 < — £b for 
A B , < 0, as is the case here. 



IV. FIDELITY 

The preceding sections demonstrate that arbitrary 
singlc-qubit rotations are possible within this set-up. The 



FIG. 3: Same as Fig.[2]except that here the rotation is about 
the 1 axis and the pulse durations are: T = 20 ps (green 
dash-dot), 50 ps (blue dashes) and 100 ps (black solid line). 
Longer times are required to perform this operation, which 
involves a transfer of population between qubit levels. Pulse 
parameters v and A were chosen to give the desired rotation 
axis and angle, A = 7T, for each value of the detuning S. 



important question then arises as to how fast these oper- 
ations can be performed. This requires that we know the 
degree to which our approximate description matches the 
actual behaviour of the system. This we quantify with 
the fidelity between the calculated operation and the re- 
sult obtained by numerical integration of the Schrodinger 
equation. The fidelity is defined as 



(20) 



where the overline represents an average over all input 
states | Vfin) , U is the predicted operation, and p out is the 
actual output density matrix. In evaluating the fidelity 
from our numerics, we have used the method described 
inRef. [l(J. 

We calculate the fidelities for the two A = it rotations 
above: one about the 3 axis and the other about the 1 
axis. In Fig. [5] we plot T as function of the laser detuning 
5 for the rotation about the 3 axis for several values of 
the pulse duration T. This figure shows that, for the 
3-axis rotation, even very short pulses (T w 5 ps) can 
be used with fidelity close to unity provided that the 
detuning S is chosen appropriately. In Fig. [3] we plot the 
same thing for the rotation about the 1-axis and observe 
that longer pulse durations are required to obtain high 
fidelities. Irregularities in T occur in both these figures 
at values of 5 = 0, As,— Y. B . These arc the regions at 
which we know from Eq. (J7]) that the perturbative theory 
breaks down. 

We now examine the dependence on the applied mag- 
netic field of the maximum obtainable fidelity for a given 
pulse-duration, . We consider a A = tt, 1-axis rotation, 
since this is the more demanding rotation, and in Fig. [4] 
plot the results for several different magnetic fields. As 
is clear, this field determines the pulse-duration required 
to obtain a given fidelity. For the parameters used here 
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FIG. 4: The maximum fidelity T for a A = n rotation about 
the 1-axis as a function of pulse duration T. Results are 
shown for magnetic fields of 2 T (green dash-dot), 4 T (blue 
dash), and B = 8 T (black solid line), with g- factors as in 
Fig. [21 Fidelity increases with both magnetic field and pulse 
duration. For B=8 T a fidelity of 95% can be obtained with 
pulse durations of ~ 50 ps. 



a field of 8 T gives a fidelity of 95% for a pulse duration 
of about 50 ps. 

The required pulse-duration is reduced if the ^-factors 
are larger. In particular, we have used here an electronic 
p-factor of \g%.\ = 0.46. This is quite small compared 
with other measurements reported for InAs QDs, which 
have yielded values of \g%\ = 0.6 @ and \g%\ w 0.9 [lj. 
Use of such dots will significantly reduce cither the mag- 
netic field or the pulse duration required to perform high- 
fidelity operations. It should also be noted that a ir- 
rotation about the 1-axis is the worst case example and 
that all other rotations can be achieved in less time. 

In any case, these time scales are much shorter than 
the trion lifetime in InAs dots which has been measured 
to be 1 ns or longer [l2|, HH. This, coupled with the fact 
that the gate operation proceeds essentially adiabatically, 
means that the effects of spontaneous emission from the 
trion has negligible effect on the gate operation. 



V. HOLE MIXING 

We now include the effects of hole- mixing in our analy- 
sis. The dominant mixing terms in the Luttinger Hamil- 
tonian [H], [l5[ mean that, rather than the 'bare' heavy- 
hole states | | , ± | ) , the actual states are better ap- 
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where I § ,=F \ 



) are light-hole states, 9 m and 4> m are 
mixing angles, and we have used the growth direction as 
the quantization axis. 

In terms of orbital and spin degrees of freedom, the 



valence-band electron states are [ljl [lj| 
3 - l±l,±±L 



2 ' 2 



2 ' - 1 - 2 



i \±l,Th)e±\ | |0,±i 



(22) 

where on the right-hand side, the first index corresponds 
to the orbital degree of freedom, 



\±l) = y/T/2(\X)±i\Y)); \0) = \Z), 



(23) 



and the second to the electron spin ±1/2. The make-up 
of these states means that, for light propagating in z- 
direction, transitions involving the light-hole components 
of \Hz) are possible and that they acquire an extra factor 
of in the matrix element. 

We are now in a position to consider the Hamiltonian 
of the system including hole-mixing, for which we use 
the same basis as for Eq. ([3]), except that now we use 
the hole-states \H±) = l/y/2(\H~) ± \H+)). If we il- 
luminate the dot with the same linear polarizations as 
before, hole-mixing means that all four transitions are 
driven by each polarization and the Hamiltonian of the 
system is far more complicated than that considered in 
the foregoing analysis. However, by adjusting the polar- 
izations of the two lasers, we can reach a situation where 
the Hamiltonian is identical to that of Eq. ([3]), but with 
renormalized parameters. The single qubit rotations then 
follow directly. 

The two laser polarizations that accomplish this are 
V = 2- 1 / 2 (cr + + e^+o-_) and and H' = 2~^ 2 (cr + - 



er_) with the phases 



•\/3cos6l m ism^e^ 
\/3cos0 m ± sin6* m eT^ 



(24) 



Note that these two polarization are non-orthogonal in 
general. Using these polarizations we obtain a Hamilto- 
nian the same form as Eq. ^ but with fly and Qh being 
replaced by 



tlx 



1 + 2cos26»„ 



= n 



3 cos 9 m 
1 



- V3e"^" 
2 cos 29 m 



smf 



H' 



3 cos B n 



y/3e i( f 



(25) 



where fiy and fin 1 are the actual Rabi frequencies of the 
two transitions. Use of these new polarizations therefore 
allows us to circumvent the effects of hole mixing and 
proceed directly as outlined in the previous sections. 



VI. GAAS DOTS 

We now describe briefly the application of the above 
theory to QDs in which the trion level does not split 
under in-plane magnetic field. This was the situation 
studied in Ref. Q, but the approach we pursue here is 
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different and ieads to an improved description of this sys- 
tem. 

With the in-plane hole g-factor equal to zero, and 
with no hole-mixing, illumination with <x + -circularly- 
polarized light propagating in the growth direction ex- 
cites a trion state \t z +) containing a spin- up heavy hole 



| fr) aligned in the growth direction. We employ two such 
<r_|_ lasers set to a detuning S from the Raman transition 
between the two electron ground states |x±). A full de- 
scription is given in Ref. [5(, where it is shown that in 
the basis {|aH-) ( \x— ), |t z +), }, the Hamiltonian for the 
system in the rotating frame is given by 



H = 
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Hi 






a iE%t—ia 
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(26) 



where Q-\(t) and are the Rabi frequencies of the two 
<T_)_-induced transitions, and E B is the Zeeman splitting 
of the electron. 

From this Hamiltonian we can derive an effective two- 
level model for the qubit sector exactly as in the previous 
sections. We obtain 



based on h^it) will typically have a fidelity of the order of 
75%. This shows that the treatment of these off-resonant 
terms is essential to the proper understanding of singlc- 
qubit rotation in QDs and that the adiabatic elimination 
approach pursued here provides just such a treatment. 
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This has the same form as Eq. © with the substitutions 



ft, 



H 



f2| and Sr, A; 



E'k. The axis and 



angle of the corresponding single-qubit rotations follow 
directly. 

Although derived in a different way. the work of Chen 
et al. Q essentially posits that the qubit is governed by 
the effective Hamiltonian 



h 2 (t) 



I^tI 2 



(28) 



This differs from the Hamiltonian hi by the absence of 
the on-diagonal terms with S ± E B in the denominator 
that describe the off-resonant effects of the qubit levels 
being driven by both lasers. These terms have a sig- 
nificant impact on the accuracy of the description with 
pulses as short as we use here. With parameters for which 
the description based on hi(t) has a fidelity of 95%, that 



VII. SUMMARY 

We have presented a theory of single-qubit rotations in 
InAs quantum dots. This theory provides a description 
of the behaviour of the system to high fidelity for realistic 
dot parameters and short laser pulses. 

Our theory shows how two Raman-detuned lasers can 
be used to obtain arbitrary single-qubit rotations, and 
that this can be accomplished with laser pulses the du- 
ration lasting a few tens of picoseconds. We have also 
shown that hole-mixing can be simply incorporated into 
this scheme through a change in laser polarizations. 

This work differs from that presented in Ref. Q not 
only in that it applies to systems in which the hole has an 
in-plane g-factor but also in the methodology. The adia- 
batic elimination approach that we have utilized here pro- 
vides better a account of the system because it includes 
the effects of the unwanted off- resonant transitions. 

This work was supported by ARO/NSA-LPS. 
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